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Abstract 

We present in this paper a calculation of the average proton-nucleus in- 
elasticity. Using an Iterative Leading Particle Model and the Glauber model 
we relate the leading particle distribution in nucleon-nucleus interactions 
with the respective one in nucleon-proton collisions. To describe the leading 
particle distribution in nucleon-proton collisions we use the Regge-Mueller 
formalism. 
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1 Introduction 



It is well known that the inelasticity is one of the most important variables to anal- 
yse cosmic ray data on hadronic cascade and on extensive air showers. A wide 
and diversified range of hadron interaction calculations in cosmic ray and accel- 
erator physics is strongly dependent on the inelasticity parameter, whose energy 
dependence is presently in question. 

Inelasticity is understood as the fraction of the available energy released for 
multiple particle production in an inelastic hadron-nucleus collision. Its value was 
estimated a long time ago from cosmic ray experiments [|I]] as being around 0.5, 
as was later confirmed at the CERN ISR [§, 

The question of the energy dependence of the inelasticity was dealt by a num- 
ber of authors through different approaches , most of them in a model depen- 
dent way. The obtained results are not consistent and the question re- 
mains unsolved. Considering the average inelasticities in hadron-proton colli- 
sions, a carefuU analysis of inclusive reaction data was done by Bellandi et al. 

. The average partial inelasticities were extracted in a model-independent way 
from experimental data of inclusive reactions initiated by pp collisions {pp — > c, 
c = n^, K^,p, p). From the results obtained for vr, K, p, p partial inelasticities it 
was also estimated the behavior of the total proton inelasticity , wich turned out 
to be rapidly increasing with energy in the high energies region. 

Bellandi et al. have already discussed this question in connection 

with the behavior of the hadronic cascade and extensive air showers in the at- 
mosphere, showing that the average proton-air inelasticity is also an increasing 
function of the energy. A model dependent analysis of the average inelasticity 
in ^ was done by means of the so called Interacting Gluon model (IGM) [|^], 
[]T0|], which includes, besides soft gluonic interactions, semi-hard QCD interac- 
tions responsible for minijet production. This model in the original version ^ 
had predicted inelasticity decreasing with energy. With the addition of the semi- 
hard component [|T^ the total inelasticity is an increasing function of the energy 

For proton-nucleus scattering, at low energy, several models for describing 
the leading particle spectrum have been proposed (Interacting Gluon model and 
Regge-Mueller formalism) [ [TTl ] , [ [T2| ] . Here, we shall work in the Iterative Leading 
Particle Model [[r3|], [|T^ and use the notation of Frichter, Gaisser and Stanev [[T5|]. 
In this model the leading particle spectrum in p + A^ A^(nucleon)+X collisions 
is built from sucessive interacions with u interacting protons of the nucleus A and 
the behaviour is controlled by a straightforward convolution equation. It should be 
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mentioned that, strictly speaking, the convolution should be 3-dimensional. Here 
we only considered the 1 -dimension approximation. 



2 Inelasticities 

In this model [ [T3| ] an iterative scheme was used to write the longitudinal dis- 
tributions for multiple scattering of nucleons with wounded nucleons. After v 
collisions the longitudinal distributions are related by means of the following 1- 
dimensional Mellin convolution integral 

Ml{x)= /'^[5jli(yK_iM^i(x/y) + 5-_i(y)(l-/3,_i)A/4%(x/y)] (1) 
Jx y 

for protons and 

M:{x) = C'^[St.M^-iM:_^{x/y) + S-_M{^-^-i)Ml_,{x/y)] (2) 
Jx y 

for neutrons. M^'"(x) are the proton and neutron distributions normalized as 

'''rfxMP'"(x) (3) 



with + nJJ = 1. The numbers express the outgoing nucleon {p and n) 
multiplicities for u wounded target nucleons. The superscripts (± ) describe in- 
teractions wich preserve and change the projectile isospin, respectively, and the 
parameters (3i, specifying the fraction of isospin preserved reactions. The S^_i{y) 
define the probability of transation of a nucleon with longitudinal momentum frac- 
tion x/y to a state with longitudinal momentum x, after v — \ collisions. For the 
probability functions Sf{y), we have for the first collision 

SHy) - ^""^^^ (4) 

with appropriate definitions of and [|T3]]. 

In this model it is assumed that we have different inelasticities upon subse- 
quent collisions. Adopting a power law form with adjustable factor for z/ > 1, we 
write 

I^dyy^^Uriy) ^ ^ 
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In order to compute the average nucleon-nucleus elasticity, < x > n-a we use 
the Glauber model [|16p. The N — A leading particle can be obtained by means of 
the relation 

oo 

Mn_a = E P^M,, (6) 

where Pj, is the probability of z/-fold collisions of the nucleon inside the nucleus, 
given by 

and 

Puib) = ^ [aZATib)rexp[-aZATib)l (8) 
where T(b) is the nuclear thickness, given in terms of the nuclear density by 

r+oo 

T{b) = / dzp{b, z) 



and normalized in the following way 

d%T{b) 



The inelastic cross section aj^ is given by the following Glauber model JT^ ] 
relation 

at^ = I d% [1 - eM-^7otAT{b)]] (9) 

As in [ [T5| ] we shall assume that the Sf{y) are the same for all interactions 
with more than one collision, u > \. Using Eqs. ([1]) and @, it is straigthforward 
to show that the serie in Eq. is an absolutely convergent serie and we can write 

(1- < X >n-a) = <K >n-a 

= -4^ f d%[l - {vexp[-{l - k)aZAT{b)] + 

+ {l-^)exp[-aZATm] (10) 

where 

k = f3,.i C dyySt{y) + (1 - dyyS-{y) (11) 



and 

n? < X >? < X >'1 < X >N 



This parameter r] defines the relationship between the nucleon elasticity in the 
first interaction and the one for the sucessive interactions of protons and neutrons 
with nucleus. The Eq. ( P^ gives a relation between the inelasticity in a nucleon- 
nucleus collision with inelasticities in nucleon-nucleon scattering. In a recent 
paper [[T^] we have used this model to analise cosmic ray data on the hadronic 
flux and we have shown that the preserved momentum fraction of the sucessive 
interactions is the same as the one for first interactions, that is k ^< x >Ar, r/ ~ 1. 
We then simply have 

(1- < X >n-a) = < K >n-a 

= -4^ I d'b [1 - exp[-(l - k)aZAT{b)]] (13) 



In this situation 1 — k =< K >jv and Eq. (|0]) gives a relationship between 
average inelasticities [ |T5| ] . 

It is clear from this relatioship that only in small aZ limit is < > 
~< K >M- In general, < K >n-a'>< K >m, the effect increasing with the 
increase of oZf If < -ft' >n^ 0, one also has < K >m-a^ 0. On the other 
hand, if < >n= ^, then < K >n-a= 1, and Eq. ( [T3| ) coincides with Eq. (|^). 

We use here the Woods-Saxon model for the nuclear distribution which 
is given by 



1 + exp[ 



(14) 

o 



where the factor {l + u^) corresponds to the Fermi parabolic distribution correc- 
tion. The parameter po is a normalization factor, 

d\p{r) = 1 (15) 

The parameters Tq, Uq and co can be derived from experimental data and we have 
To = 0. 976^4^/^ fm, Oq = 0.546 fm and the parameter cu is given by 

Lo = -0.25839 if A < 40 



00 = if A > 40 

In order to calculate < K > n^a we use for < K >j\f the values calculated 
by means of the Regge-Mueller formalism [ [T2| ] and as input for aZ we have used 
the UA4/2 parametrization for the energy dependence [|19]]. In the Fig. (1) we 
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show the resuks of this calculations for the following nuclei: C, Al, Cu, Ag, Pb 
and air (A=14.5). In this figure we also show recent emulsion chamber data for 
p-Pb, <K >= 0.84 ± 016 ^] and for p-C, < K >= 0.65 ± 0.08 [0. 

In the Fig. (2), we compare the calculated < K >p-air with results from some 
models used in Monte Carlo simulation [|22|]; the Kopeliovich et al. [ ^3| ] (KNP) 
QCD multiple Pomeron exchanges model; the Dual Parton model with sea-quark 
interaction of Capella et al. ; the statistical model of Fowler et al. [ |25| ] and 
with calculated values derived from cosmic ray data by Bellandi et al. . We 
note that the calculated < K >p-air in was done assuming for the T(6) 
nuclear thickness the Durand and Pi model p7|], which gives small values for the 
average inelasticity. In the Fig. (2) we also show the average inelasticity values as 
calculated by means of this model. 



3 Conclusions 

We have here calculated the average proton-nucleus ineslaticity in the Glauber 
framework, relating the leading particle distributions in nucleon-nucleus interac- 
tions with the respective one in nucleon-proton collisions. We have compared our 
results with recent emulsion chamber data for p — Pb [20] and for p — C [ ]2T| ] at 
piab = 1.20 X 10^ GeV/c. At least in the experimental errors limit our calcula- 
tion is in agreement with these experimental data. We have also calculated the 
average p — air inelasticity in a wide range of energy. In order to describe the 
nuclear thickness we have used two models: the Woods-Saxon model [[TSp and 



the Durand and Pi model [^. The average proton-air inelasticities calculated by 
means of the Woods-Saxon model are larger than the ones calculated by using of 
the Durand and Pi model. 

One remark should be stressed. The calculated < K >p_ajr values derived 



from cosmic ray data pq ] were obtained assuming an approximation for the lead- 
ing particle distribution in proton-air collisions. Therefore, it is model dependent. 
The discrepancies between the values of the < K >p-air at low ^/s are conse- 
quence of the fact that two different sets of experimental data were used: nucle- 
onic flux and hadronic flux at sea level ( for discussions see p^). Finally, we note 
that the calculated < K >p-air with the Woods-Saxon model shows a behavior 
with ^/s wich goes between the values calculated by means of the QCD multiple 
Pomeron exchanges model (KNP) and that one calculated by means of the Dual 
Parton model. 
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Figure 2: The < K >^ as a function of ^/s in GeV. The experimental data 
from [E^. Dashed line from [|25|]. Full line from [p4]]. Dot-dashed line from [[23|]. 



Dotted line from Eq.([T3|) with Woods -Saxon model [fl^]. Long dashed line from 
Eq. (|13j) with Durand-Pi model 
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